Molecular Phylogenetics and Evolution

Animal mitochondrial rRNAs are commonly viewed as structurally reduced, yet sponge mt SSU rRNAs range from compact to highly expanded structures. Using nine conserved structural anchors, we compared 216 taxonomically resolved records from four classes and 22 orders, including 16 freshwater Spongillida and 200 marine sponges. Twelve homologous hypervariable substructures were coded as structural types, and their ordered combinations as composite types. We identified 38 structural types and 62 composite types across molecules ranging from 853 to 2,019 nt. Hexactinellida and freshwater Spongillida were each uniform for a distinct composite type but differed markedly in overall structure: hexactinellid mt SSU rRNAs were compact, whereas those of Spongillida were long and contained four to five candidate insertion regions. These results show that a conserved scaffold can accommodate extensive lineage-associated structural variation and provide a practical framework for comparing highly divergent mitochondrial rRNAs.

Determining the correct donor species/lineages of grass-to-grass lateral gene transfer (LGT) is vital for deducing specific donor features that could help inform the mechanism of transfer. This requires a dataset spanning a broad range of species to achieve the phylogenetic resolution necessary for precise donor inference. As grass-to-grass LGT often involves the transfer of multi-gene DNA fragments, they can contain additional sequences that allow for accurate orthologous comparisons, such as nuclear DNA of plastid origin (NUPTs). Here we systematically scan for NUPTs in the genomes of four Alloteropsis semialata accessions, whose LGTs have previously been characterised. Using the abundant Panicoideae chloroplast sequences, we reconstruct NUPT phylogenies and infer two lateral acquisitions: one from Paniceae/Digitaria and another from Andropogoneae/Eremochloa adjacent to a previously identified LGT. We then assembled and included an additional 12 Eremochloa chloroplast genomes in the analysis and showed the likely donor was Eremochloa attenuata. Subsequent short-read mapping from E. attenuata to the nuclear region flanking this NUPT showed consistent coverage across the region, including the previously identified LGT, supporting co-transfer. Overall this study highlights the potential for NUPTs to better identify the donors of grass-to-grass LGT.

Protists play important roles in food chains and symbioses in soil and aquatic environments, displaying an enormous morphological and functional diversity. While most commonly found protist species are well known to science, our global-scale environmental DNA survey across soil, water, and sediments reveals dozens of novel, phylum-level phylogenetic lineages that remain to be characterized for basic morphology and function. A vast majority of these undescribed taxa occur in marine water and sediments, but some are common in soil. Most of these novel taxa have distinct substrate and habitat preferences and biogeographic patterns. To accord these lineages scientific agency and enable unambiguous scientific communication, we propose formal names for 150 species to phylum-level taxa from 25 deep lineages based on eDNA and rRNA gene long-read sequence information.

Sicyoidochytrium minutum DNA virus (SmDNAV) was isolated several years ago from a protist host of family Thraustochytriaceae of the class Labyrinthulomycetes. This virus shared little similarity to other viruses in gene content and protein sequences, albeit seemingly belonging to the phylum Nucleocytoviricota. By extensive searches in genomic and metagenomic sequence databases, we identified numerous long contigs related to the SmDNAV genome and analyzed proteins shared by these putative viruses. Phylogenetic analyses place these viruses within the class Megaviricetes, outside of all established orders, and as a sister group to the clade combining families Mamonoviridae and Manesviridae. Homologs of SmDNAV proteins were found in association (either integrated or co-sequenced) with other Labyrinthulomycetes and Rhodophyta protists from diverse marine and freshwater environments. Consequently, we propose SmDNAV as the prototype member of a new order, provisionally named Ariadnavirales, within class Megaviricetes, phylum Nucleocytoviricota. Members of Ariadnavirales have lost most of the genes encoding components of the replication and transcription systems that are otherwise conserved in nucleocytoviricots, suggestive of transition to genome replication and expression dependent on the host nucleus.

Background Trypanosoma cruzi, the causative agent of Chagas disease, is subdivided into distinct genetic groups known as Discrete Typing Units (DTUs), each with distinct genetic traits that influence epidemiology and transmission dynamics. Several triatomine species serve as potential vectors of T. cruzi in the United States. However, despite the growing number of Chagas disease cases in the country, little is known about the genetic diversity and population structure of T. cruzi in natural vector populations. Methodology/Principal Findings We applied a multilocus metabarcoding approach to improve DTU resolution and characterize the genetic diversity and structure of T. cruzi in triatomines collected across five states of the southern United States. Five single-copy nuclear markers and one mitochondrial marker were amplified and processed by high-throughput sequencing to assess genetic diversity. We recovered 35 nuclear and 15 mitochondrial haplotypes from 70 infected specimens. Overall, genetic diversity was low ({pi} < 0.01 at all nuclear loci), with DTUs TcI and the North American lineage of TcIV detected, TcI being the most prevalent. Geographic structuring was particularly evident in TcI strains, which exhibited a distinctive haplotype profile in Florida populations, potentially linked to the recently revalidated vector species Triatoma ambigua. Mitochondrial introgression from TcIV into TcI suggests inter-DTU genetic exchange in these populations. Multiple haplotypes within individual insects detected across single-copy nuclear markers, support multiclonal infection as common feature of T. cruzi in natural vectors. Conclusions/Significance These findings provide new insights into the genetic landscape and evolution of T. cruzi in the United States. Evolutionary connectivity through mitochondrial introgression and frequent multiclonality highlights the importance of deep sequencing approaches for resolving T. cruzi genetic diversity, with direct implications for understanding for transmission dynamics, disease monitoring and control.

Hydrozoa is a group of relatively simple animals with a well-developed nervous system. The nervous system in all hydrozoan medusae is highly conserved and includes outer and inner nerve rings at the bell margin, a neural network in the manubrium, and radial neural pathways that connect them. However, one element of the nervous system shows substantial variability among species: the subumbrella neural network. We examined the structure of the nervous and muscular systems in the subumbrella of 14 species of hydrozoan medusae. The main conclusion of this study is that the distribution of neural networks in the subumbrella strongly correlates with the distribution of smooth radial muscles. This correlation suggests that smooth radial muscles are the primary target of the subumbrella nervous system. Most species in the order Anthoathecata show a trend toward secondary loss of the neural networks and radial smooth muscle fibers in the subumbrella region, concentrating neural elements and smooth muscles only in the radial pathways along the radial canals. By contrast, all studied species in the order Leptothecata have neural networks in the subumbrella area, as well as radial smooth muscle fibers spread throughout the entire subumbrella region. The correlation between radial smooth muscles and the nervous system is also observed in the radial pathways along the radial canals. All species with thick bundles of smooth radial muscles along the radial canals have clearly defined, dense neural pathways running along or even embedded within the smooth muscle bundles.

Aegilops peregrina is a wild allotetraploid wheat wild relative and an important source of genetic diversity for stress tolerance and agronomic traits. Here, we report a subgenome-resolved, chromosome-scale reference genome assembly of a drought tolerant and stem rust resistant Ae. peregrina accession PI 604178 generated using PacBio HiFi and Hi-C sequencing. The 10.13 Gb assembly contains 98.81% of sequence anchored to 14 pseudomolecules representing the seven S and seven U chromosomes, with contig and scaffold N50 values of 25.84 and 746.48 Mb, respectively. The assembly achieved a consensus quality value of 74.61, 97.83% k-mers completeness, and 99.9% BUSCO completeness. LTR Assembly Index values of 20.43 and 18.79 for the S and U subgenomes, respectively, further supported high continuity across repeat-rich regions. Repetitive elements comprise 85.93% of chromosome-anchored assembly. We annotated 59,910 high-confidence protein-coding genes, with comparable gene representation across the two subgenomes. This reference genome provides a high-quality genomic framework for comparative analyses, characterization of important loci regulating agronomic and resilience related traits, and sequence-guided exploitation of Ae. peregrina allelic diversity for wheat improvement.

Biological nitrogen fixation converts dinitrogen gas into ammonia, supplying new bioavailable nitrogen to marine ecosystems, but the evolutionary processes shaping its distribution among heterotrophic bacteria remain unresolved. Thalassolituus, a genus within the family Oceanospirillaceae (order Oceanospirillales), is best known for hydrocarbon degradation, yet nitrogen fixation has been confirmed in only one cultured isolate. We analyzed 74 quality-filtered genomes assigned to Thalassolituus within a broader dataset of 421 Oceanospirillaceae genomes to reconstruct the distribution and evolutionary history of the minimal nifHDKENB gene set. Twenty-five genomes encoded complete or near-complete nif loci and occurred in four well-supported clades interspersed with genomes lacking the pathway. Statistical topology tests rejected the species-tree topology for concatenated NifHDK and NifHDKENB protein alignments, and eleven recombination events across nif loci were supported by at least four detection methods. The core nifHDK gene order remained broadly conserved, but accessory neighborhoods differed among clades, and structural nifHDK genes showed stronger codon adaptation than biosynthesis nifENB genes. Clade 2 combined species-gene tree congruence, conserved gene neighborhoods, and comparatively high nifH codon adaptation, whereas Clades 1 and 4 showed greater phylogenetic discordance, more recombination, and weaker codon adaptation. These results support a reticulate history in Thalassolituus, in which lateral acquisition introduced nitrogen fixation into distinct lineages, vertical inheritance preserved it within some clades, and homologous recombination continued to reshape nif loci. These processes help explain why nitrogen fixation is unevenly distributed among closely related marine heterotrophic bacteria.

The study of evolution between species (phylogenetics) and the study of evolution within a species (population genetics) are highly related, as the same biological mechanisms are fundamental to both fields. Although both have been studied for a long time, their joint study in a unified setting has been prevented by the different time scales they consider and the different data types they employ. A similar discrepancy holds for their whole-genome specializations, comparative genomics and pangenomics. Two active areas in these fields are genome rearrangement studies and graphical pangenomics, respectively. Since the emergence of graphical pangenomics, these have existed as separate fields, despite observations that central data structures representing genomic variants in both fields are highly similar. While there exists a wealth of theoretical results for various rearrangement models in comparative genomics, the application to pangenomic data is hampered by the limitations of rearrangement problem formulations. On the practical side, pangenomes typically contain too many individual genomes for classical problems, such as the often NP-hard parsimony problems, to be solved, or for all-vs-all comparisons using rearrangement distances to be performed. On the theoretical side, some assumptions in the formulation of rearrangement problems, such as the assumption of an underlying tree, are inadequate for many pangenomes. In this work, we propose the Complete Ancestral Reconstruction for Pangenomes (CARP) problem, which overcomes these limitations while retaining intuitive relationships to both classical rearrangement problems and pangenome graphs.

Motivation: The Interactive Tree of Life (iTOL) is widely used to display and annotate phylogenetic trees, but managing its format-sensitive annotation files impede reproducible high-throughput analyses. Among the maintained Python packages and versions evaluated, none combined template generation, taxonomic monophyly assessment and iTOL batch operations. Results: PyiTOL validates inputs, generates 31 iTOL template schemas (22 accepted by the live batch uploader), performs LCA-based monophyly classification with nested-monophyly detection, sampling-completeness states and polyphyletic subgroup decomposition, plus API upload and session replay. On a topology-constructed benchmark, all calls matched prespecified labels for 4,389 groups; on a 700-genome tree, binary mono/non-mono calls agreed with ETE4 for 409 genera; 17,294 GTDB R232 genera were processed in about 17 s. Availability and Implementation: PyiTOL 1.0.3 (Python [&ge;]3.10; Linux, macOS and Windows) is MIT-licensed at https://github.com/ZengZichao/PyiTOL and archived with test data at Zenodo (https://doi.org/10.5281/zenodo.22106806).

Extrafloral nectaries (EFNs) are a widespread plant defense mutualism trait and are highly convergent, appearing in hundreds of plant lineages worldwide. Here we investigate a report of possible EFNs in Erythranthe angulosa, a recently-described California wildflower. We integrate field observations, insect bioassays, an induction experiment, and microscopy to test for signatures of EFN function, finding no evidence that the distinctive axillary swellings produced by E. angulosa function as EFNs. We also uncovered two distinct morphs at the type locality of E. angulosa that diverge in the number of axillary swellings produced, as well as other shoot architecture traits such as stem thickness, leaf size, and branch number. Although the axillary swellings appear to not function as EFNs, they remain a compelling morphological variant within the yellow monkeyflowers that may perform storage or another unknown function.

Genomic signatures can provide key insight into the evolutionary history and remaining adaptive potential of threatened populations. As demographic decline erodes both diversity and the processes maintaining it, understanding how remaining variation is distributed becomes increasingly important for conserving species like the staghorn coral, Acropora cervicornis, a foundational but critically endangered Caribbean reef-builder. We analyzed 46 high-coverage A. cervicornis genomes from 10 locations across the tropical western Atlantic to evaluate neutral and adaptive structure, genomic diversity, demographic history, inbreeding, and connectivity, and generated a regional haplotype reference panel for future genomic monitoring. Genome-wide analyses recovered recurring regional substructure, but differentiation was modest and partly explained by isolation-by-distance and spatial variation in effective migration. Subpopulations had similar levels of genomic diversity, shared demographic history, and limited evidence of local adaptation. These patterns support interpreting sampled Caribbean populations as a single evolutionarily significant unit (ESU) containing multiple regional management units (MUs), rather than as deeply divergent evolutionary lineages. Despite substantial retained variation and low current inbreeding, estimated contemporary effective population size was small, suggesting an increased vulnerability to the effects of drift as demographic collapse continues, especially if structure is reinforced by isolated management. Together, our findings emphasize the urgent need for interventions that preserve and enhance genomic diversity, including risk-managed assisted gene flow. Supported by the haplotype reference panel developed here, these strategies will require coordinated efforts across regional entities to conserve and restore A. cervicornis as a jointly managed, single ESU.

Although the biosynthetic pathways of benzylisoquinoline alkaloids (BIAs) have been extensively investigated in several plant species, their transcriptional regulatory mechanisms remain only partially understood. Jasmonate (JA)-responsive group IX APETALA2/Ethylene Responsive Factor (AP2/ERF) transcription factors (TFs) are well-known regulators of specialized plant metabolism, including the biosynthesis of various alkaloids. However, their specific roles in BIA biosynthesis remain largely elusive. Here, we isolated five novel group IX AP2/ERF TFs, designated Benzylisoquinoline alkaloid Jasmonate-responsive AP2/ERF (BJE1-5), from Coptis japonica. Phylogenetic analysis revealed that Benzylisoquinoline alkaloid Jasmonate-responsive AP2/ERF (BJE) proteins belong to subclades distinct from group IXa, which contains well-known AP2/ERF TFs involved in alkaloid biosynthesis. Transient expression analyses in C. japonica protoplasts demonstrated that certain BJEs, particularly CjBJE3 and CjBJE5, positively regulated BIA biosynthetic genes through a mutual regulatory network among BJE members. Moreover, CjBJE3 expression was regulated by CjbHLH1, a unique-type basic helix-loop-helix (bHLH) TF specific to BIA-producing plants. Furthermore, heterologous expression of CjBJE3 and CjBJE5 in cultured Eschscholzia californica cells significantly enhanced the overall BIA production, particularly by increasing end-product benzophenanthridine BIAs, highlighting several uncharacterized biosynthetic genes clustered in the genome. Our findings suggest that BIA-producing species have developed a specific regulatory network comprised of CjbHLH1 and BJE TFs, providing valuable clues for identifying novel biosynthetic enzymes.

High-phenolic plant genera present a major technical limitation in genomic research. Standard extraction approaches that perform reliably across diverse flora often perform poorly when applied to recalcitrant taxa, producing low DNA yield and integrity incompatible with sequencing requirements. The genus Anigozanthos (Kangaroo paws) from the family Haemodoraceae exemplifies this problem. We identified key physicochemical factors governing extraction failure in this genus and resolved them through targeted modifications to lysis chemistry and contaminant management. The resulting protocol achieved a near threefold improvement in DNA purity, substantially reducing contaminant carry over and consistently yielded high-integrity, long DNA fragments (DIN > 7) across a diverse sample set spanning cultivated and wild material across four diverse genera of Haemodoraceae. We also tested a straightforward purity assessment framework that can be implemented in any standard molecular laboratory, enabling rapid pre-submission quality assessment without the need for specialised equipment. Together these advances open a practical path to genomic characterisation of Anigozanthos that establishes a transferable model for genomic research across Australia ' s chemically complex native flora.

Candida auris (currently Candidozyma auris) is an emerging fungal pathogen responsible for dramatic global increase in invasive candidiasis with high mortality. Most concerning, C. auris has a high propensity to colonize patients and persist and develop multidrug resistance to main classes of antifungals. In this study, we investigated the genetic and phenotypic diversity and resistance mechanisms of C. auris clinical isolates recovered from hospitalized infected patients. A total of 53 isolates from 38 unique patients were recovered from various clinical sources and evaluated for susceptibility to routine antifungal drugs. Whole genome sequencing (WGS) and single nucleotide polymorphism (SNP) analysis were performed to generate a phylogenetic network to infer population structure and identify mutations associated with drug resistance development. Isolates were also phenotypically evaluated for ability to form biofilms and aggregate, and cell wall adhesins gene expression studies were performed to provide mechanistic insights into C. auris phenotypic plasticity. Except for one clade III isolate, all isolates belonged to clade I and all were resistant to fluconazole with incidence of resistance to amphotericin B, echinocandins or both. Non-synonymous SNPs were found in genes associated with antifungal resistance including ERG11, TAC1B, CDR1 and FKS1. Phenotypically, isolates varied in their ability to form biofilm and aggregate which correlated with expression of the Scf1 and Als4112 cell wall adhesins genes highlighting C. auris phenotypic plasticity in circulating clinical strains. These findings underscore the growing clinical threat posed by C. auris and reinforce the need for optimized surveillance and treatment strategies for controlling its spread.

Objective: This study examined the prevalence of psychotic experiences (PE) and how early onset and persistence of PE contribute to risk and severity of mental disorders in adolescents at familial high-risk of schizophrenia (FHR-SZ) or bipolar disorder (FHR-BP) and adolescents from a population-based control group (PBC). Methods: This is the second follow-up of a nationwide cohort study including 522 children at FHR-SZ (N=202), FHR-BP (N=120), and PBC (N=200). Participants were assessed at ages 7, 11, and 15 using a semi-structured interview to evaluate PE and mental disorders. Results: At age 15, adolescents at FHR-SZ reported more PE than PBC over the past six months (current) and the past four years, while adolescents at FHR-BP only reported more current PE. PE reported at two or three timepoints (persistent PE) predicted any Axis I disorder in mid-adolescence, corresponding to three- (OR 2.9, 95% CI [1.5-5.7]) and 21-fold (OR 21.4, 95% CI [2.8-162.3]) increased risks, respectively. Persistent PE also predicted multimorbidity, with three- (OR 2.8, 95% CI [1.0-7.6]) and four-fold (OR 4.1, 95% CI [1.2-14.1]) increased risks, respectively. This was after adjustment for sex, early mental disorders, and familial risk. Conclusions: This study demonstrates a strong link between persistent PE and mid-adolescence mental disorders. Our findings emphasize PE as important risk markers for mental disorders during mid-adolescence and highlight the importance of monitoring children with PE before age 7 who develop persistent symptoms.

Next-generation sequencing (NGS) for the detection of somatic variants has become the method of choice in a variety of molecular oncology fields and in the clinic. Its use ranges from sequencing entire tumor genomes and transcriptomes to targeted clinical diagnostic gene panels. The NYU Langone Genome PACT (Profiling of Actionable Cancer Targets, LG-PACT) assay is a qualitative in vitro diagnostic test that uses targeted next generation sequencing (NGS) of formalin-fixed paraffin-embedded (FFPE) tumor tissue matched with normal specimens from patients to detect gene alterations in a targeted panel covering 606 genes and the TERT promoter. Indications for testing are cancer (solid tumors and hematological malignancies) where a mutational profile from multiple genes would be informative for disease stratification, prognosis, or treatment options including targeted therapies and eligibility for clinical trials. The test is intended to provide information on somatic mutations including point mutations, small insertions/deletions (indels), and copy number aberrations for diagnostic and treatment decisions. LG-PACT is a United States Food and Drug Administration (FDA) cleared diagnostic test (510K: K202304). The clinical interpretation of sequencing data of molecular tumor markers from NGS encompasses automated variant calling tools with human interpretation. This final mostly manual review of data step is intensive, involving highly trained scientists, encompassing literature review, interpretation and clinical tier classification by pathologists, who then provide a complete molecular diagnostic report to the treating oncologists. We provide analysis of 1339 clinical genomic profiles from 31 different cancers and their subtypes, comprising of central nervous system (CNS) 792 (59%) cases (incl. meningioma, glioma and glioblastoma), with 267 (20%) cases predominantly of lung, pancreatic and colorectal and 280 of others (21%). Here, we present the technical challenges of validating an NGS oncological diagnostic targeted assay for clinical grade accuracy and sensitivity for patient care. We show how copy number alterations provide a more comprehensive description of the tumors genomic profile. We then outline the utility of targeted panel sequencing based on certified pathologist selection of reportable variants for our current patient cohort. Where analysis of variant detection has led to 49.4% (661/1339) of our clinical tumor samples containing mutations in known therapy targeted genes, 35.6% (477/1339) with mutation detected in other genes, and 15% (201/1339) cases being negative.

Importance: Patients are increasingly learning Alzheimers disease (AD) genetic and biomarker results through electronic health portals. Evaluation of alternative scalable delivery models for return of AD risk information is needed to best support patient understanding and psychological well-being. Objective: To determine whether a patient-centered digital platform is comparable to clinician-mediated telehealth sessions for returning APOE and plasma pTau-217 results on outcomes of knowledge and psychological well-being. Design: The Evaluation of Self-Mediated Alternatives for Risk Testing Education and Return of Results (eSMARTER) study was a noninferiority trial of a patient-centered digital platform compared to clinician-mediated disclosure of APOE genotype and optional pTau-217 disclosure. Setting: Decentralized, fully remote trial enrolled participants in the contiguous United States (U.S.) between October 2024 and February 2025, with follow-up completed in November 2025. Participants: Eligible participants were aged 60-80 and had previously undergone APOE genotyping (without disclosure) via the GeneMatch program, passed psychological screening, had internet access, and were English-speaking. Interventions: Participants were randomized, 2:1, to the eSMARTER digital platform or clinician-mediated disclosure of APOE genotype. Following the 6-month post-APOE assessment, participants were offered optional pTau-217 disclosure via the same randomized modality. Main Outcomes and Measures: Primary outcomes at 1-7 days following APOE disclosure included changes in anxiety, disease-specific distress, and AD-related knowledge within a priori non-inferiority margins. Results: 674 persons (mean [SD] age 68 [4.7] years; 451 [67%] female; mean [SD] telephone MoCA=19 [2]) were eligible and provided demographic information. 651 participants were randomized to clinician-mediated (n=216) or digital disclosure (n=435) and completed APOE disclosure (66 [10%] APOE4 homozygotes, 377 [58%] heterozygotes, 208 [32%] non-carriers). 604 participants completed the study; 500 completed optional pTau-217 disclosure. Baseline characteristics were balanced across groups. At 1-7 days following APOE disclosure, scores on AD-related knowledge, PROMIS Anxiety, and disease-specific distress measures met non-inferiority. Conclusions and Relevance: Disclosure of APOE genotype by the eSMARTER digital platform is non-inferior to clinician-mediated telehealth disclosure. No significant between group differences were found following disclosure of pTau-217 results. Together, these results suggest that this digital platform may provide an evidence-based scalable approach for returning AD genetic and biomarker results.

Abstract: Objective: Post-traumatic epilepsy (PTE) is a common sequela of traumatic brain injury (TBI). Research indicates that individuals with PTE tend to experience greater cognitive difficulties compared to those with TBI alone. However, it is plausible that a distinct cognitive profile exists that distinguishes between TBI cases with and without PTE. We aimed to identify longitudinal changes in cognitive measures among TBI patients to better assess the changes associated with developing PTE. Setting: Outpatient. Participants: Prospective subjects who had suffered TBI within 6 months post-injury (TBI-6M, n=32), retrospective subjects with pre-existing PTE diagnoses (PTE, n=20), and healthy control subjects (HC, n=41). Design: We examined cognitive performance for TBI patients within 6 months post-injury, then again within 12 months (TBI-12M, n=26), and within 18-months (TBI-18M, n=25), and compared this with cognitive performance among HC and PTE. Main Measures: Cognitive tests administered yielded 15 test components for analysis. We utilized linear mixed effects modeling to examine cohort-level differences cognitive function. Results: 11/15 tests showed a significant performance deficit in the PTE subjects compared to HC. TBI-6M was not significantly different from the PTE subjects; with time, 9/15 tests showed some degree of recovery in TBI subjects. Tests for information processing speed/working memory and executive function showed strong recovery (TBI-6M vs. TBI-18M, SDMT written: p<0.0001, SDMT oral and COWAT: p<0.001). Tests for visual attention/working memory also showed a smaller but significant recovery (TBI-18M vs. PTE, p<0.05). By contrast, tests for verbal memory [HVLT-R Delayed Recall] showed chronic impairment in TBI (TBI-18M vs HC, p<0.0001). TBI subjects generally trend towards recovery in cognitive performance post-TBI. Conclusions: Information processing speed/working memory are strong indicators for TBI recovery, while auditory learning/memory shows chronic impairment. The stagnation of recovery in cognitive domains typically characterized by robust recovery may correlate with an elevated risk of developing PTE.

The distribution of fat-soluble compounds (including antioxidants) in follicular fluid (FF) remains sparsely documented in relation to in vitro fertilization (IVF) outcomes and existing studies have reported diverging associations. This study aimed to describe plasma and FF concentrations of fat-soluble micronutrients in women undergoing IVF and to analyze their adjusted associations with ovarian function, embryo development and pregnancy outcomes. In 2021-2022, plasma and FF samples were collected from 82 women (first IVF cycle) at oocyte puncture, along with lifestyle data covering the three preceding months. Eleven compounds (two tocopherols, three xanthophylls, five carotenes and retinol) were quantified. All compounds were detected in both compartments (lowest in FF) except phytoene, undetectable in FF. Plasma and FF -tocopherol concentrations were positively associated with plasma estradiol levels before oocyte puncture (both p<0.01) while FF -carotene and lycopene were inversely associated with plasma progesterone concentrations (p=0.01 and 0.02, respectively). Plasma phytofluene and phytoene were positively associated with mature oocyte rate (p=0.03 and p=0.01, respectively), while FF retinol was negatively associated (p=0.03). Carotenes, tocopherols and retinol were inversely associated with later IVF outcomes: fertilization rate (p<0.001 for plasma g-tocopherol, 0.02 for FF retinol), top-quality embryo (p=0.02 for plasma phytofluene), biochemical pregnancy at day 7 post-embryo transfer (p=0.05 for plasma -tocopherol, 0.02 for plasma -carotene), clinical pregnancy (p=0.03 for plasma -tocopherol, 0.01 for plasma phytoene) and live birth (p=0.04 for plasma -tocopherol, 0.02 for plasma phytoene). Plasma and FF g-tocopherol were positively associated with embryo fragmentation (both p<0.05). Finally, among xanthophylls, only plasma {beta}-cryptoxanthin was positively associated with plasma progesterone concentrations (p=0.02). Our findings of heterogeneous associations between tocopherols, carotenes, retinol and IVF outcomes across the stages of IVF suggest a beneficial effect limited to early outcomes and support a complex and context-dependent role of these compounds in female reproduction. This manuscript has been submitted to PlosOne on August 19, 2026.